Breathers in T'T-symmetric optical couplers 
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We show that the parity-time (VT) symmetric coupled optical waveguides with gain and loss 
support localised oscillatory structures similar to the breathers of the classical 4> 4 model. The 
power carried by the PT-breather oscillates periodically, switching back and forth between the 
waveguides, so that the gain and loss are compensated on the average. The breathers are found to 
coexist with solitons and be prevalent in the products of the soliton collisions. We demonstrate that 
the evolution of the small-amplitude breather's envelope is governed by a system of two coupled 
nonlinear Schrodinger equations, and employ this Hamiltonian system to show that the small- 
amplitude PT-breathers are stable. 
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I. INTRODUCTION 

Light propagation in PT-symmetric optical systems 
with balanced gain and loss has been under intense 
scrutiny in the past few years. The concept has its 
roots in quantum mechanics where a VT symmetric non- 
Hermitian Hamiltonian may have an entirely real spec- 
trum of eigenvalues 0, 0. In optics, the VT symme- 
try can be achieved by an appropriate modulation of the 
complex refractive index [3j-l5|. 

The symmetric optical systems should display a va- 
riety of unusual and often counter-intuitive phenomena 
including an unconventional beam refraction [g, 0] , Bragg 
scattering || nonreciprocal Bloch oscillations ficj . 
symmetry-breaking transitions [ill [l2j . a loss- induced 
optical transparency 13 1, the conical diffraction [l3 |. a 
new type of Fano resonance [l5[ , chaos [l6| , and nonlocal- 
ity manifested in the nontrivial effect of the boundaries 
[l7j . Recently, optical PT-symmetric couplers [i~2l. Il3T] 
and lattices [18| have been realised experimentally. 

Nonlinear effects in PT-symmetric systems are of par- 
ticular interest for the fundamental and applied science. 
They offer potential for an efficient control of li ght, in- 
cluding the all-optical low-threshold switching |19l - [21j 
and unidirectional invisibility (20j . In addition, nonlin- 
earity can compensate the diffraction of stationary light 
beams and dispersion of light pulses allowing the forma- 
tion of spatial and temporal solitons. 

There has already been a large number of studies of 
optical solitons in PT-symmetric systems. Solitons in 
complex one-dimensional potentials were analyzed on the 
basis of the nonlinear Schrodinger equation [22| - l32j . The 
two-dimensional symmetric potentials were dealt with in 
Refs. [H H, H. The authors of [H]-[37| classified soli- 
tons in the planar PT-symmetric couplers, whose geom- 
etry is intermediate between one- and two-dimensional 
lattices. 

The PT-symmetric solitons considered in the above 
publications represented stationary self-localised modes. 



The solitons arise due to the exact compensation of the 
gain and loss at each moment of time. A more general 
type of localised objects was identified in [37[ where the 
unstable solitons were observed to seed spatially-localised 
temporally-periodic states. (In the context of planar sta- 
tionary waveguides, these are interpreted as the transver- 
sally localised structures with profiles oscillating along 
the waveguide.) These objects resemble breathers in con- 
servative systems (such as the </> 4 and sine- Gordon equa- 
tion) [38j ; hence they were referred to simply as breathers 

M- 

In this paper the VT breathers are studied in more 
detail. First, we derive the amplitude equations for the 
oscillatory solutions in the planar PT-symmetric nonlin- 
ear optical coupler (equations for the envelopes of the 
oscillatory wavepackets). The amplitude equations turn 
out to be Hamiltonian — despite the fact that the orig- 
inal system includes both gain and loss. These Hamil- 
tonian equations are then used to show that the (zero- 
velocity) PT breathers form two-parameter families with 
variable amplitude, localisation width, and contrast of 
power density oscillations. We also employ these equa- 
tions to establish the stability of the breathers with small 
amplitude. Finally, the planar PT-symmetric coupler is 
simulated numerically. Results of our numerical simula- 
tions demonstrate that the breathers are generic objects 
which are commonly formed as a result of the soliton 
collisions. 

The outline of the paper is as follows. In Sec. [Til we in- 
troduce the mathematical model, and in the subsequent 
section, derive equations for the slowly- varying envelopes 
of its oscillatory solutions. Section IIVI uses these ampli- 
tude equations to classify the PT-symmetric breather 
states. The stability of the small-amplitude breathers 
is established in section [V] In Sec. I VII we describe the 
formation of breathers in the soliton-soliton collisions. 
Finally, Sec. IVIII summarises results of this study. 
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nor their sum are conserved. The total power satisfies 



FIG. 1. (Color online) A schematic representation of VT- 
symmetric coupled waveguides with gain (red) and loss (blue 
waveguide), (a) Two planar waveguides carrying stationary 
light beams. Here t and x indicate the longitudinal and 
transversal spatial coordinate, respectively, (b) A pair of one- 
dimensional waveguides where light pulses undergo temporal 
evolution as they travel along the x axis. 



II. MODEL 

The PT-symmetric coupler, i.e., a pair of coupled 
waveguides with power gain in one waveguide and op- 
tical loss of equal rate in the other, has been studied the- 
oretically 0, i, IM IH EH and experimentally [H [lj] . 
Optical systems that include the PT-symmetric coupler 
as a structural element [HI, Il7l . lioll a nd s yste ms consist- 
ing of arrays of such couplers [Jill H M, M, EJ |H| 
have also been discussed in literature. 

Following [35|, [13] we analyze the diffraction of optical 
beams propagating in a planar VT coupler, in media 
with the Kerr-type nonlinearity. The amplitudes of the 
active and passive modes in this setting satisfy a system 
of two coupled nonlinear Schrodinger equations, 



iut + u xx + 2\u\ 2 u = —v + i'yu, 
ilk + v xx + 2|u| 2 t> = — u — ijv. 



(2.1) 



Here t is the (spatial) coordinate in the propagation di- 
rection and x is the transversal coordinate. The coeffi- 
cient 7 > is the amplification rate for the waveguide 
with gain and, at the same time, the damping rate for 
the waveguide with loss. This planar coupler is schemat- 
ically shown in Fig. Ufa), ft i s fitting to note here that 
the system (|2 . 1 1) emerges as the continuum limit of the 
chain of VT couplers considered in [34j . 

The same VT symmetric system (|2.ip can describe 
the propagation of optical pulses (rather than station- 
ary light beams) [37| . This alternative interpretation of 
Eqs. (|2.1[) arises if t and x stand for the time and distance 
in the frame of reference travelling along with the pulse. 
This is the arrangement illustrated by Fig. Hfb). 

The system (|2.1j) is not conservative. Neither the in- 
dividual powers associated with the two modes, 



V u = 



~dx, 



(2.2) 



j t (V u +V v ) = 2 1 {Vu-V v ), 



(2.3) 



which implies that it remains constant only on solutions 
which have V u = V v for all times [37]. 



III. WEAKLY NONLINEAR AMPLITUDE 
EQUATIONS 

We start our analysis by transforming Eqs. (|2.1[) to a 
system with a diagonal linear part. Assuming 7 < 1 and 
defining 



2uj 



where 



v = arcsm7, 
Eqs. (|2.1j) are taken to 

idt + a xx - cj a 4 
+4ie- lt> ja 2 b* 
ibt + b xx + uj b - 



2uj 
uj{) = cos 9, 

2(|a| 2 + 2|6|> 
I- 2e 2i0 a*b 2 = 0, 
2(2|a| 2 + |6|> 



(3.1) 



-4ie i0 ^a*b 2 + 2e~ 2i V&* = 0. 



(3.2) 



The system Q3.2p has two simple reductions or, equiva- 
lently, two invariant manifolds. Letting b = 0, Eqs. (|3.2j) 
reduce to a scalar nonlinear Schrodinger equation 



ia t 



Lo^a + 2\a\ 2 a = 0, 



(3.3) 



while letting a = yields a scalar Schrodinger equation 
with the opposite sign of the frequency term: 



ib t + b xx + uj b + 2\b\ 2 b = 0. 



(3.4) 

Both (|3.3p and (|3.4[) have soliton solutions and hence the 
system (|3.2j) admits two types of 'simple' solitons: one 
with 6 = and the other one with a — 0. These low- 
and high-frequency solitons have been analysed before 
[34I l35l. [37| . Here, our aim is to construct more general 
solutions with both components nonzero. 

To this end, we note that when a and b are so small that 
the nonlinear part in (|3.2p can be neglected, the result- 
ing linear system has a family of spatially homogeneous 
stationary- wave solutions: a = Aoe~ luJot , b = Boe lu>ot . 
To search for the nonlinear counterparts of these, we con- 
sider a long-wavelength small-amplitude configuration: 

a(x, t) = e 1/2 A{X, t), b(x, t) = e 1/2 B(X, t), (3.5) 

where X = e 1 / 2 ^ and e a small parameter (e > 0). The 
0(1) fields A and B satisfy 



iA t + eA xx - lo q A + 2e(\A\ 2 + 2\B\ 2 )A 

+4ie- l8 e~fA 2 B* 
iB t + eB X x + u B - 

-Aie w ejA*B 2 - 



+ 2e zw eA*B z = 0, 
2e{2\A\ 2 + \B\ 2 )B 
- 2e- 2l \A 2 B* = 0. 



(3.6) 
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Solutions of the system (13.6[) can be sought for as ex- 
pansions in powers of e: 



A = A + eAx 



B = B Q + eBi 



(3.7) 



We also assume that the coefficients A n and B n depend 
on a hierarchy of 'slow times' and 'zoomed out' spa- 
tial coordinates: A n — A n (To,Tx, ...;Xo,Xx, ...), B n = 
B n (To,Ti, ...;Xq,Xi, ...), where 



T n = e n t, X n = e n X, n = 0,1,2. 



(3.8) 



In the limit e — > the scaled time and space variables 
decouple, and can be treated as independent. In what 
follows, we adopt a shorthand notation 

D n = d/dT n , d n = d/dX n . 

Note that the parameter e is not pegged to any scale 
of the original model (|2.1|) . (|3.2p . Therefore we expect it 
to be absorbable in the parameters of solutions that we 
will end up with. 

Substituting the expansions (|3 . T[) in (|3.6p . we equate 
coefficients of like powers of e. The order e° gives 

(iD - uj )A = 0, 
(iD +Lu a )Ba = 0, 



whence 



with 



A = e lT p, B = e lT q, 



(3.9) 



t = wo T . 

The coefficients p and q are functions of all variables ex- 
cept Tq. 

The order e 1 produces 

{Wo - wo)^ - -[WxAo + d 2 A Q + 2(\A \ 2 + 2\B Q \ 2 )A 

+4ie- te 1 A 2 B* + 2e 2t6 B 2 A*}, 

{iD + uj )B 1 = -[iD x B + d 2 B + 2{\B \ 2 + 2\A \ 2 )B 

-4ie ie ^B 2 A* + 2e- 2l6 A 2 B*}. 

(3.10) 



To eliminate the secular terms, we impose 

iDx P + d 2 p + 2(\p\ 2 + 2\q\ 2 )p = 0, 
lD iq + d 2 q + 2(\q\ 2 + 2\p\ 2 )q = 0. 



(3.11) 



The remaining terms in the right-hand sides of (|3.10[) 
involve the third harmonics only; hence we get, for Ai 
and Bi, 



a 1 / 2i9 2 * 3ir a ■ -i8 2 * -3ir\ 

A\ = - — (e q p e — 4ie jp q e ) 

Bx = -— Uie^ iq 2 p*e^ +e- 2i9 p 2 q*e- 3ir ) 
2uj " 



(3.12) 



Proceeding to the order e 2 , and setting the correspond- 
ing secular terms to zero, we obtain 



iD 2 p + 2d d lP + —(\q\ 2 - 2\p\ 2 )\q\ 2 p = 0, 
iD 2 q + 2d dxq + —(2\q\ 2 - \p\ 2 )\p\ 2 q = 0, 

UJQ 



(3.13) 



where we have substituted for Ax and Bx from (|3.12[) . 

According to Eqs. p.lip . the variations in the ampli- 
tudes p and q become noticeable only over long periods of 
time, At ~ e _1 . Eqs. (|3.13|) govern the evolution of these 
amplitudes over even longer time intervals, At ~ (T 2 . It 
is convenient to combine Eqs. (|3.1ip and p,13p into a sys- 
tem that takes care of the evolution on both slow scales. 
To this end, we add Eqs. ([3TTj) to Eqs. (l3~T3l multiplied by 
e and define T = et. Since the amplitudes do not depend 
on T , the chain rule gives d/dT = D x + eD 2 + e 2 D 3 + .... 
Thus, to within C(e 2 ), we have Dxp + eD 2 p = pt and 
Dxq + (D 2 q = qx, and so the resulting pair of equations 
can be written as 

IPT+PXX + 2(\p\ 2 + 2\q\ 2 )p + —(\q\ 2 - 2\p\ 2 )\q\ 2 p = 0, 
iq T + qxx + 2(\q\ 2 + 2\p\ 2 )q + -(2\q\ 2 - \p\ 2 )\p\ 2 q = 0. 

UJQ 

(3.14) 

(Here e > 0). This is a hamiltonian system, with the 
Hamilton functional 



H = J [\px\ 



\qx\ 



111 



A\pq\ 2 ) 
>■)} dX. 



The amplitude equations p. 141) describe the evolution 
of the slowly changing envelope of a small-amplitude, 
weakly localised packet of waves with the carrier fre- 
quency wo- Over time intervals e _1 < At < e~ 2 , equa- 
tions (|3.14j) are equivalent to the original system (|2.ip . 
This remarkable equivalence of a dissipative and conser- 
vative system, holding for a particular but fairly broad 
class of trajectories, is attributable to the 'PT-symmetry 
of the former. 

Setting e = 0, the system p.!4p becomes 



IPT+PXX + 2(\p\ 2 +2\q\ 2 )p = 0, 
iqr + qxx +2(\q\ 2 +2\p\ 2 )q = 0. 



(3.15) 



This vector nonlinear Schrodinger equation has been ex- 
tensively studied in literature [43j-|55[. On the other 
hand, the system (|3.14p with £ / does not seem to 
have been discussed before. 

Note that both Eq. p.l4p and the "curtailed" system 
p,15p govern the small- amplitude breathers only, with 
u, v ~ e 1 / 2 . However Eq. p,14p has an advantage over 
Eq. p,15p in that the former system has a longer range 
of validity. While Eq. p,15p ceases to be valid for times 
exceeding e _1 , Eq. p,14p remains accurate for times as 
long as e~ 2 . 
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Another reason for the evaluation of the second or- 
der corrections in the perturbation expansion, is related 
to the conservativity of the amplitude equations (|3.14j) 
and p,15[) . Once the first-order amplitude equations are 
found to be given by a hamiltonian system [the system 
p. lip ], the question arises whether this property is spe- 
cific to the first-order evolution only. The fact that the 
second-order dynamics are also governed by a hamilto- 
nian system, suggests then that the conservativity is an 
inherent property of the small-amplitude oscillations. We 
conjecture that this property is valid to all orders in the 
perturbation theory (and may only be violated by terms 
that lie beyond all orders). 



IV. BREATHER SOLUTIONS 

One simple solution of Eqs. (|3.14p is 

p = e <^+¥*)^ sec h[y/I(X- VT)}, q = 0. (4.1) 
The other one is given by 

p = 0, q = e i ( vT+ ^ x )^sech[Vv(X - WT)]. (4.2) 

These two solutions of (13.141) will be referred to as de- 
generate solitons. The parameters > 0, v > 0, V and 
W can be chosen arbitrarily. Here and v give the am- 
plitudes of the degenerate solitons, and V, W are their 
velocities. 

The degenerate soliton solutions of Eq. (|3.14p corre- 
spond to the solitons of the scalar reductions (|3 . 3[) and 
(I3.4p of the original system (|2.ip . The degenerate soliton 
(14. ip corresponds to the low- frequency soliton of (|2.1I) . 
and the solution (|4.2p to its high-frequency counterpart 
[MH! S3- The vector of the power densities {\u\ 2 , \v\ 2 } 
associated with each of these solutions describes a pulse 
travelling, without oscillations, at the velocity v = e 1 / 2 ^ 
and w — e 1 / 2 ^, respectively. 

Our main interest is in solutions of the system (I3.14p 
which have both components nonzero. Thanks to the 
Galilian invariance of (|3 . 14[) . it is sufficient to con- 
sider separable solutions corresponding to nonpropagat- 
ing waves: 

p = e^ T P(X), q = e tuT Q(X). (4.3) 
The spatial parts P and Q satisfy 

P" - fiP + 2{\P\ 2 + 2\Q\ 2 )P + — (|Q| 4 - 2\PQ\ 2 )P = 0, 

Q" -vQ + 2(\Q\ 2 + 2\P\ 2 )Q + — (2\PQ\ 2 - |P| 4 )Q = 0, 

wo 

(4.4) 

where we use the notation " = d 2 /dX 2 . 

Localised solutions of the stationary equations (I4.4p 
give rise to oscillatory, breather-like, configurations in 




FIG. 2. (Color online) Numerical evolution of the initial con- 
dition in the form of the expansion (I4.12p - (|4.14p with t — 0. 
(In this simulation, 7 = 0.5 and e = 0.1.) Shown are \u\ 2 (a), 
I?; 1 2 (b), and powers carried by the two components of the 
breather (c). In (c), Vu is depicted by broken red and V v by 
dotted blue line. Also shown is the total power V u +V V (solid 
line). The simulation continued until times much longer than 
e -2 = 100, without any visible change in the amplitude or 
period of the breather. 
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the original model (|2.ip : 

= ,1/2 



v(x, t) 



+ 0(e 3 / 2 ), 



u(x, t) = e 1 ' 2 [Qe l " 2t + Pe^'J + 0{e A/z ), 

c l/2 \Q e i(L0 2 t+9) _ p e -i{u x t+ff) 

where P = P{e 1 / 2 x), Q = Q{t x l 2 x), and 

LOi = UJq — 6(1, Ll>2 = UJq + (V 

The corresponding |u| 2 and |i>| 2 are 



e(|P 

2 



QP*e 



* i(uil+LU2)t 



+ C.C. 



e (|P| 2 + |Q| 2 )- e QP*e 



where c.c. stands for the complex conjugate of the im- 
mediately preceding term and we neg lected the 0(e 2 )- 
corrections. These quantities show temporal oscillations 
with the frequency lo\ + u>2 = 2 cos 9 + e(v — /i). 

In this paper, we confine ourselves to the simplest 
choice of v — \i. (A brief comment on a more general sit- 
uation with v 7^ /i is in the Appendix IA1) An additional 
simplification is attained by restricting to real solutions. 
For real P and Q equations (|4.4|) reduce to 



P" - /xP + 2P 3 + 4Q 2 P + — (Q 2 - 2P 2 )Q 2 P = 0, 
Q" - fiQ + 2Q 3 + 4P 2 Q + — (2Q 2 - P 2 )P 2 Q = 0. 

UJQ 

(4.5) 



where we have defined £ = fi x ^ 2 X. 

The operator in the left-hand side of (|4.8p has a zero 
eigenvalue, with the associated eigenfunction being odd. 
If we wish to construct a solution with definite parity 
(i.e. an even solution), we should take 5 = 0. On the 
other hand, the operator in the left-hand side of (14.9[) is 
positive definite, hence invertible. As a result, Eq. (|4.9[) 
has an exponentially decaying solution: 



1 u 3/2 

V = =- — (6sech£ + sech 3 £)- 

51V3 wo 



Taken together with 5 = 0, this implies 



1 u?l 2 

n = -Pi = ,= - — (6sech£ + sech 3 £). (4.10) 



102^3 wo 



Returning to the original variables u and v we note 
that, as expected, the parameters e and [i enter the solu- 
tion only in combination e/i. Without loss of generality, 
we can set one of these to 1, e.g. /i = 1. 



Equations (|4.3p . with P and Q expanded as in (|4.7|) . 
and P n , Q n as in (|4.6I) . (|4.10l) provide solutions to the 
amplitude equations (13.141) : 



When e = 0, the system (|4.5p has an explicit solution 

Po(A) = Q Q (X) = ^| sech(^A). (4.6) 

The terms proportional to e in (I4.5[) are regular pertur- 
bations, i.e., the perturbed solution satisfying the bound- 
ary conditions P(X),Q(X) — > as \X\ — > oo exists for 
all sufficiently small e. To show this, we expand P and 
Q in powers of e, 



P = P + ePi + e 2 P 2 



Q = Qo + eQi + e 2 Q 2 



(4.7) 

and substitute the expansions in (|4.5p . Letting S = P\ + 
Qi and T> = Qi — Pi, the order e gives 

(-d 2 /d£ 2 + l-6sech 2 £)S = 0, (4.8) 



p = — — sechA 
v/3 

o = — = sechA 
V3 



1 - 



102a; 



102w 



(6 + scch 2 A) + 0(e 2 ) 



(6 + scch 2 A) + 0(e 2 ) 



(4.11) 



Since both p and q are nonzero in (I4.1ip . we will be re- 
ferring to these solutions as two-component solitons. 

Feeding Eqs. (I4TT1) in ([33]), J37|, flU]), (l3~T2|) gives 



a = e 1 / 2 [A + eAi + 0(e 2 )] , 
6 = e 1 / 2 [Bo + eBx + 0(e 2 )] , 



(4.12) 



with 
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g-i(wo— e)t 

A = 7 = sech(e 1/2 x) 



B 



V3 

e i(ui +e)t 



V3 



Ai 
Bx = 



6 y/3 0Jq 



sech(e 1 ^ 2 x 
sech^e 1 / 2 



(6+ scch 2 (e 1/2 x)) + 0(e 2 ) 
102.0 (6+-chV^)) + 0(e 2 ) 



102w 
e 



3 i(3u t+2fl) 



0(6), 



6-v/3wo 



sech^e 1 / 2 ^) 4i7e i(3^ot+e) + e -i(3 W0 *+2fl) +0(e) _ 



(4.13) 



Equations (14. 121) - (14. 131) . taken together with the conver- 
sion formulas 



u(x, t) = a + b, v(x, t) = e l6 b — e 



(4.14) 



yield solutions of the original equation (|2.1[) . 

To test the accuracy of the asymptotic solution (|4.12[) - 
(|4.14j) . we simulated equations (|2.ip with the initial con- 
ditions in the form (|i7l2|) - (|4TTi|) with t = 0. [In these 
initial conditions, we neglected the 0(e 2 ) terms in Aq, Bq 
and the 0(e) terms in A\,B\\ The resulting oscillatory 
configuration is plotted in Fig. [2J The fundamental har- 
monic in the frequency spectrum of \u\ 2 and |u| 2 was in- 
deed found to be very close to 2wo, the double frequency 
of the asymptotic solution. 

As we mentioned in section UH the system (|2.ip may 
be thought of as a continuum limit of a chain of coupled 
PT-symmetric dimers. The power in each dimer can per- 
form a periodic oscillation [23, [2l[ , with an amplitude- 
dependent period. The breather is an oscillation involv- 
ing the entire chain. Although the amplitude of oscil- 
lation varies along the chain, the coupling synchronises 
individual dimers so that the breather has a single base 
frequency. Accordingly, the power integrals (|2.2p associ- 
ated with the two modes show a perfectly periodic be- 
haviour [Fig. Htc)]. 

The total power V u + V v is not a constant of motion 
but is periodic and therefore, conserved on average. 



STABILITY 



The amplitude equations (|3.14[) may be used to study 
the dynamics of the solitons and breathers of the original 
system (|2.1[) over times up to t ~ e~ 2 . In particular, 
Eqs. (|3.14p may be used to study the stability of these 
objects. 

Consider a stationary solution (|4.3p of the system 
(|3.14p . This can be one of the two degenerate soli- 
tons (|4.ip and (|4.2p — or the nondegenerate soliton 
(|S7T] ) .p~5 J) , (|130]) corresponding to the breather of the 
original system (I2.ip . We consider the simplest situa- 
tion where /i = v; in this case we may set, without loss 
of generality, fi = v = 1. Linearising Eqs. (I3.14p about 
the stationary solution and assuming perturbations of the 



form 

5 P (X,T) = e lT [{(X, T) + ig(X, T)] 
5q(X,T)=e lT [y(X,T)+i z (X,T)] 

where f , g, y and z are real, gives 

£if+V(X)y=-g T , £ g = f T , 
Miy + V(X)i = -z T , M z = y T - 

Here we have introduced the operators 



(5.1) 



C = 


-d 2 /dX 2 h 


- 1 


- 2P 2 


-4Q 2 H 


- — (2P 2 

UJQ 


-Q 2 )Q 2 


A = 


-d 2 /dX 2 4 


- 1 


- 6P 2 


-4Q 2 H 


-— (6P 2 

UJQ 


-Q 2 )Q 2 


M a = 


~d 2 /dX 2 - 


- 1 


-4P 2 


-2Q 2 - 


- — (P 2 - 
UJQ 


2Q 2 )P 2 


Mi = 


-d 2 /dx 2 - 


- 1 


- 4P 2 


-6Q 2 H 


- — (P 2 - 


6Q 2 )P 2 



UJQ 



and a coefficient function 



4e 



V(X) = -8PQ + — (P 2 - Q 2 )PQ. 

UJQ 

For separable solutions of the form 

f(A, T) = Re [e XT f(X)] , g(A, T) = Re [e XT g(X)] , 
y(X, T) = Re [e XT y(X)] , z(X, T) = Re [e XT z{X)\ , 

with complex f,g,y,z, and A, Eq. (|5.ip reduces to an 
eigenvalue problem: 



where 



X 



Ho 
-Hi 



(5.2) 



(5.3) 



is a 4 x 4 matrix with blocks given by 
Ho = 

and y, z are two-component vectors 



Co 
Mo 



_( Ci V(X)\ 
V(X) Mi J' 



V 






z = 











(5.4) 
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A. Stability of the high- and low-frequency solitons 

Consider, first, the degenerate soliton (I4.ip and let the 
velocity V = 0. [This degenerate soliton with Q = de- 
scribes the amplitude of the low-frequency soliton of the 
original PT-symmetric equations (|2.1[) .] In this case, the 
operators Cq and L\ reduce to L and L l7 respectively, 
where 

L Q = -d 2 /dX 2 + 1-2 scch 2 A, (5.5) 
Li = -d 2 /dX 2 + 1-6 scch 2 A, (5.6) 

while Mo and M.\ acquire a common form which we 
denote Li: 



Li = -- 



dX 2 



1 - 4 sech 2 A + — sech 4 A. (5.7) 



tv 



Since Q = implies V(A) = 0, the eigenvalue problem 
(I5.2[) acquires a block-diagonal form: 



Log 

L i y ; 



A/, Li/ = 
— Az, i 1 2 



Ay. 



(5.8) 
(5.9) 



Eq. (|5.8p is the linearised eigenvalue problem for the 
scalar cubic nonlinear Schrodingcr equation, a well re- 
searched integrable system. It has no discrete eigenval- 
ues except the four-fold zero eigenvalue. Its continuous 
spectrum occupies the imaginary axis. 

On the other hand, Eq. (|5.9[) gives 

L\v = -A 2 y. 

2 

This implies that A = iuj, where w is an eigenvalue of the 
hermitian operator Li. Since all such eigenvalues are 
real, all A's are pure imaginary and hence the degenerate 
soliton is stable. 

When e = 0, the operator Li has two discrete eigen- 
values, u! a and cub, given by 

w a = a-3« -1.438, u b = 3a - 4 « 0.685, (5.10) 



with a = (VT7— l)/2. The corresponding eigenfunctions 
are ip a = sech" A and ipb = sech Q_1 Atanh A, respec- 
tively. The eigenvalues uj a and Wf, persist when e deviates 
from zero (but remains finitely small). It is only when 
e grows above a certain finite value that Wf, and then ui a 
immerse in the continuous spectrum. Accordingly, for e 
below a finite threshold, the degenerate soliton (|4.1j) has 
two modes of internal oscillation. (For e = 0, this fact 
has been established in [Hj].) 

The degenerate soliton (|4.2[) corresponds to the high- 
frequency soliton of the original equations (12.11) . The 
linearisation about this degenerate soliton leads to the 
same eigenvalue problem (|5.9I) . with the same operator 
(I5.7[) . where one just needs to replace e — > — e. This 
observation establishes the stability of the soliton (|4.2I) . 
As long as e remains below a finite threshold, the operator 



Li with e — ► — e has two discrete eigenvalues; hence the 
degenerate soliton (|4.2[) has two internal modes. 

The fact that the degenerate solitons of the amplitude 
equations p.!4|) are stable implies that both the low- 
and high-frequency solitons of the PT-symmetric system 
(I2.1[) are stable for sufficiently small e. This conclusion 
is in agreement with the analysis of the low- and high- 
frequency soliton performed directly on the equations 
(|2.1[) . Namely, the high-frequency soliton was shown to 
be stable when its amplitude a lies below a finite thresh- 
old a c , a c = (|) 1/2 (1 - 7 2 ) 1 / 4 [H H3- On the other 
hand, the low-frequency soliton has an unstable eigen- 
value irrespectively of the amplitude but its real par t 
is exponentially small when the amplitude is small [37j . 
This instability constitutes an effect that lies beyond all 
orders in e"; it cannot be captured by the amplitude 
equations (|3.14|) . The unstable perturbations take an 
exponentially long time to grow in this case; hence the 
small-amplitude low-frequency soliton will not reveal any 
instability when studied over time intervals t ~ e~ n . 

The frequencies of the internal modes of the low- and 
high-frequency soliton solutions of Eqs. ()2.1[) were also 
computed in [37l]. These coincide with the frequencies 
(|5.10[) computed using the amplitude equations (|3.14[) . 



B. Stability and spectrum of breather: t ~ e 



Turning to the two-component soliton (|4.11j) . it is con- 
venient to consider the soliton of the "curtailed" sys- 
tem (|3 . 1 5(1 first. The stability of the soliton of the sys- 
tem (|3.15p will imply the stability of the breather of the 
original PT-symmetric system (|2.1|) over time intervals 
t ~ (where e 1 / 2 is the amplitude of the breather). 



The two-component soliton of the system p,15j) is 
given by Eqs. (14. lip with e = 0: 



n/3 



i tT sechA, q 



JT 



sechA. 



(5.11) 



Depending on the context, this symmetric solution was 
referred to as the linearly polarised [48| or equally mixed 
51] soliton. Note that setting e = in Eqs. (14. lip does 
not mean that we are considering breathers of zero am- 
plitude. The nonzero parameter e remains present in the 
corresponding breather solution (I4.12p . (14.131) . (|4.14p ; in 
particular the amplitude of the breather remains equal 
to e 1 / 2 + 0. 

The stability of the soliton (|5.11l) was proved by the 
construction of a Lyapounov functional [46J. With an 
eye to addressing the situation of general e, we reconsider 
the stability of this solution here — using the eigenvalue 
analysis. 



When e = 0, the eigenvalue problem (15.21) can be cast 
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in the block-diagonal form 





L 


L 



-Li 




-L 



+ 



Ci 
Pi 

P2 



Ci 

c 2 

Pi 

P2 



where the operators Lq and L\ are as in ([5 



L + = — 



tfA 5 



1 



- secli A. 
3 



(5.12) 
(5.13) 
and 
(5.14) 



In fact it is not difficult to show that the operator 
(|5.13[) does not have any discrete eigenvalues at all - 
neither real nor imaginary. (See the Appendix iBl) The 
implication is that when e = 0, the two-component soli- 
ton of the vector nonlinear Schrodinger does not have in- 
ternal modes. (This fact has been previously established 
by numerical means [53j.) With regard to the breather 
of the 'PT-symmetric system (|2.ip , this implies that the 
small-amplitude breather cannot have any modulating 
frequencies of order e in its spectrum. This is the second 
conclusion of this subsection. 



The components of the column vectors in (|5.12j) - (|5.13j) 
are the sums and differences of the components of the 
vectors in (|S~4|) : Ci = z + .9, C2 = V + f, pi = z - g, 

P2 = y - /■ 

Eq. (|5.12j) arose in the previous section [see Eq. (|5.8j) ]. 
It is the linearised eigenvalue problem for the scalar cu- 
bic nonlinear Schrodinger equation. As discussed there, 
the matrix-differential operator (|5 . 12[) does not have any 
discrete eigenvalues except the four zeros. Therefore 
Eq. (|5.12|) can be safely disregarded and we can focus 

on Eq. dnunj)- 

In order to transform Eq. (|5 . 1 3[) to a form more 
amenable to analysis, we note that the only discrete 
eigenvalue of the operator (15.14)) is f3 + 1/3, where (3 = 
v/11/12 — 1/2 > 0. (It is associated with the nodeless 
eigenfunction ip — sech^X.) Hence the operator L + is 
positive definite and admits an inverse. This observation 
allows us to write the vector equation (|5.13[) as a gener- 
alised eigenvalue problem for a pair of scalar operators, 



L0P1 



(5.15) 



In (I5.15p . Lq is a symmetric operator, and sym- 
metric and positive definite. All eigenvalues (—A 2 ) of the 
problem (|5 . 1 5[) are real and the corresponding eigenfunc- 
tions can also be chosen real. The lowest eigenvalue, 
can be found as the minimum of the Rayleigh quotient 



\2 



A = min 



(pi,L pi) 
(pi,L+ l pi)' 



(5.16) 



Here (, ) stands for the scalar product in the 
space of square integrable real functions: (4>, -0) = 
^(X)^{X)dX. 

The lowest eigenvalue of the Schrodinger operator Lq 
is zero; it is associated with the nodeless eigenfunction 
z^(X) = sechAT. Therefore the Rayleigh quotient in 
(|5.16p cannot take negative values and its minimum is 
exactly zero: — Aq = 0. This means that the matrix- 
differential operator in the left-hand side of (|5.13l) does 
not have any nonzero real eigenvalues A and so the soliton 
(15. lip of the vector nonlinear Schrodinger (|3.15l) is stable. 

This is the main conclusion of this subsection. It 
implies that the small-amplitude breather of the PT- 
symmetric system (I2.ip is stable over time intervals t ~ 



C. Stability of the breather: t ~ e 



To extend the breather stability result to times of order 
e~ 2 , we need to consider the system (|3.14p with e ^ 0. We 
should demonstrate that its solution (|4.1ip does not have 
unstable eigenvalues with Re A of order e a , < a < 1, in 
its spectrum. 

We begin the stability analysis of this solution with the 
identification of symmetries of the system (|3.14l) . These 
will provide information on zero eigenvalues of the oper- 
ator (|Q)1 . 

Besides the translation and Galilean invariance, the 
system l|3.14p is symmetric with respect to the U(l) x 
U(l) transformations of the form p —> pe 1 ^, q — > qe lx , 
where 4>,X = const. In addition, fi and v can be chosen 
arbitrarily in the stationary system (|4.4I) . Thus each so- 
lution of the form (|4.3p is a member of a six-parameter 
continuous family and therefore, the eigenvalue problem 
(|5.2p has six zero eigenvalues. 

The corresponding eigenvectors and generalised eigen- 
vectors of the matrix J? can be found explicitly. First, 
we observe that 



— Ho 



= 0, 



(5.17) 



and -Hi(P x ,Qx) T = 0; hence (P,0,0,0) T , (0,Q,0,0) T , 
and (0,0, Px,Qx) T are the U(l) and translational 
eigenvectors, respectively. One can also check that 
Ui{P^Qn) T = ~(P0) T , Hi(P„Q,) T = -(0,Q) T and 



Uow = (Px,Qx) T , where 



±X(P,Q) T . These 



define the generalised eigenvectors: (0, 0, P^, Qfj) , 
(0,0,P„g,) T , and -iX(P,Q,0,0) T . 

All nonzero eigenvalues A of the matrix J? can be found 
from the solution of the eigenvalue problem for a 2 x 2 
matrix: 



HqHi 



(5.18) 



Using (|5.17l) one can readily check that the eigenvectors 
of H0H1 corresponding to — A 2 7^ satisfy 



f(X)P(X)dX = / y(X)Q(X)dX = 0. 
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These orthogonality constraints define a subspace of the 
space of square integrable vector-functions. On this sub- 
space, the operator "Ho admits an inverse and (|5.18j) can 
be written as 

The components P(X) and Q(X) of the solution (14. 7[) 
remain positive for all X as long as e remains small. This 
means that zero remains the lowest eigenvalue of the op- 
erators Co and AAq — the operators whose null eigenvec- 
tors are given by P and Q. Therefore, the operator Hq 1 
remains positive definite (and symmetric) — while the 
operator 'Hi is symmetric. Eq. (I5.19[) implies then that 
all eigenvalues (—A 2 ) are real, so that all A are either real 
or pure imaginary. 

As e grows from zero, the six eigenvalues of the ma- 
trix R remain at the origin. New discrete eigenvalues 
can only arise by bifurcating from the continuous spec- 
trum which fills the imaginary axis of A outside the gap 
(— i, Once an eigenvalue has detached from the con- 
tinuum, it can move along the imaginary axis toward the 
origin. However the eigenvalue could only reach the ori- 
gin as e exceeded a finite threshold. Therefore, the two- 
component soliton will remain stable as long as e remains 
small. 

Concerning the breather solution of the system (|2.ip , 
the implication of this result is that the PT-symmetric 
breather is stable on the timescale t < e~ 2 . (That is, the 
breather's lifetime is no shorter than e~ 2 ). 

VI. BREATHER PRODUCTION IN SOLITON 
COLLISIONS 

Breathers are known not to be exceptional or isolated 
occurrences in the 'PT-symmetric planar coupler. In par- 
ticular, they form as a result of the soliton instability 
[37l l56j . In this section we argue that breathers are even 
more common than solitons themselves: a collision of a 
high- and a low-frequency soliton produces two or more 
breathers, and a collision of two breathers also results in 
one or more of these oscillatory objects. 

We use Eqs. (|2.1I) to simulate the evolution of the ini- 
tial condition in the form of two solitons of equal ampli- 
tudes, moving toward each other with equal velocities: 

u(x, 0) = a + b, v(x, 0) = e t0 b — e~ l9 a, 

a = e i ^ x+x °) v fjlsech[^(x + xo)], (6.1) 
b = e-^-^y/JIsechl^/flix - x )]. 

Taking xo > 0, the low-frequency soliton is initially on 
the left and the high-frequency one is on the right; the 
initial velocities are v and —v, respectively. [Note that in 
(|6.ip , the same symbol v denotes the velocity of the soli- 
ton and the second component of the vector field, v(x, t); 
this slight abuse of notation should cause no confusion.] 



The initial distance between the solitons is assumed to 
be much larger than their widths: 2^/Jlxq 3> 1. 

The high-frequency soliton is stable if fi < | yl — 7 2 
[35l [37j ■ The low- frequency soliton is unstable for all /1 
but when the amplitude is small, its instability growth 
rate is exponentially small in fi [37| . Therefore when 
the solitons' amplitudes are sufficiently small, the low- 
frequency soliton will not manifest instability in the run- 
up to the collision. The two small-amplitude solitons can 
be considered as two stable entities. 

The collision of the low-frequency and the high- 
frequency solitons in the VT symmetric system (|2.ip cor- 
responds to the collision of degenerate solitons (|4.ip - (l4.2p 
governed by the amplitude equations (|3.14p . In the par- 
ticular case e = 0, such collisions were studied by Tan 
and Yang (55[ (see also [5l|). Depending on the solitons' 
initial velocities, the colliding degenerate solitons were 
recorded to pass through each other or bounce back. The 
solitons emerging from the collision would no longer be 
degenerate; instead, they would have both p and q com- 
ponents nonzero [5l|, [55[ . Translated in the language of 
the PT-system (|2.ip . this means that the collision of the 
small-amplitude PT solitons should typically result in 
the emergence of two breathers. 

This is indeed the scenario that we have observed in our 
numerical simulations of Eqs. (|2.ip . We have detected 
the formation of two breathers in collisions of small- and 
moderate-amplitude solitons. A typical evolution is de- 
picted in Fig. [3] 

An interesting feature of the degenerate-soliton colli- 
sions recorded by Tan and Yang [55| . was that the re- 
duction of the collision velocity would not result in the 
decrease of the velocities of the solitons after collision. In 
agreement with this amplitude-equation effect, our sim- 
ulations of the collision of VT solitons with initial veloc- 
ities v — > have produced breathers diverging at finite 
speeds (see e.g. Fig. |4ja,b)). 

Another inelastic effect detected in the curtailed ampli- 
tude equation (I3.15p . pertained to the initial velocities in 
the range 0.1 < V < 0.3. For these V, the collision of two 
degenerate solitons was seen to result in the production 
of a stationary small-amplitude soliton, in addition to 
the two transmitted or reflected ones [55| . A similar phe- 
nomenon accompanies the collision of the low- and high- 
frequency small-amplitude solitons in our PT-symmetric 
system (|2.ip . Namely, the initial condition (16. ip with v 
in the range 0.1/W 2 < v < OS^, 1 ^ 2 and small /i gives rise 
to three breathers. Two of these move apart while the 
third, small-amplitude, breather is left behind near the 
origin. We have observed this effect even for not-very 
small soliton amplitudes, Fig. SJa). 

As the amplitudes of the colliding solitons are in- 
creased, the curtailed equation (|3.15p ceases to be appli- 
cable. The collision of larger-amplitude solitons is now 
accompanied by intense radiation, while the oscillations 
of the emerging breathers acquire a low-frequency mod- 
ulation [Fig. Hfb)]. As the amplitudes exceed a certain 
threshold, the collision results in a blowup of one of the 
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FIG. 3. (Color online) The collision of the low- (initially on 
the left) and the high-frequency soliton (initially on the right). 
As the solitons approach each other, they develop the beat- 
frequency oscillations of growing amplitude. The localised ob- 
jects emerging from the collision remain oscillatory despite the 
growing separation distance — these are a pair of breathers. 
The breathers are weakly radiating; also note the emission of 
a rapid small-amplitude breather at the moment of collision. 
In this simulation, 7 = 0.5, y/JI = 0.3, v — 0.4, and xq = 16. 



fragments. 

One more range of parameter values where the equa- 
tion p. 151) does not furnish any accurate description of 
the dynamics, pertains to large v. As v is increased, we 
observe the growth of the transient amplitude of one of 
the emerging breathers — a kind of a rogue wave ap- 
pearing just after the collision [Fig. Etc)]. Eventually, 
this rogue wave seeds the blow-up of the breather. 

It is worth emphasising here that the creation of 
breathers is characteristic only for the collision of two 
solitons of different types (that is, collision of the low- 
with the high-frequency soliton). The scattering of two 
like solitons, e.g. two high-frequency solitons, is purely 
elastic — for the simple reason that the initial condition 
and the resulting solution belong to the same invariant 
manifold a — 0. The constraint a = defines a reduction 
to a completely integrable equation [Eq. (I3.4|) ]. hence the 
elasticity of collisions. 

The ubiquity of the breathers stems from the fact that 
they are not confined to the a = or b = manifolds. 
They represent trajectories evolving out of generic ini- 
tial conditions which do not belong to either of the two 
reductions. 

Finally, we touch upon the collision of two breathers. 
The outcome of this collision can be predicted on the 
basis of the amplitude equation p.!4[) . Indeed, the scat- 
tering of two generic solitons in a Hamiltonian system 




FIG. 4. (Color online) The collision of solitons with moderate 
and large amplitudes, small and large initial velocities. In (a), 
7 = 0.6, y/JL — 0.3, and the initial velocity v — 0.075 lies in 
the interval (O.l/i 1 ^ 2 , 0.3pi 1//2 ). Note a small-amplitude non- 
propagating breather left behind while two large-amplitude 
fragments shoot out of the collision. The panel (b) shows the 
collision of solitons with larger amplitudes. Here 7 = 0.5, 
yfp. = 0.5, and v — 0.125. The panel (c) corresponds to large 
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FIG. 5. (Color online) The collision of two breathers. Both 
breathers are taken in the form (I4.12|l . (|4.13|) . (|4.14|) . with the 
amplitudes i/e = 0.3, and Galilei-boosted with the velocities 
v — ±0.5. The panels (a) and (b) are different in the initial 
phase of the breathers. In both plots, 7 = 0.3. 



typically produces two solitons of lower energy, or their 
bound state. Consistently with these expectations, the 
numerical simulations of Eqs. (|2.1j) demonstrate the pro- 
duction of one or two breathers (Fig. [5]). 

VII. CONCLUDING REMARKS 

Stationary solitons in the VT symmetric planar cou- 
pler are known to be sustained due to the exact offsetting 
of the power gained in the active waveguide by the power 
lost in its passive counterpart [34, 35, 37]. In this paper, 
we have described another realisation of the gain-loss bal- 



ance, which is provided by the breathers. In the breather 
case, the total power is conserved not at every moment 
in time, but only over a period of oscillation. 

Results of our study can be summarised as follows. 

1. We have derived a system of amplitude equations 
[Eqs. (|3.14[) ] governing the envelope of the breather. For 
times t < e~ 2 , where e 1 / 2 gives the scale of the ampli- 
tude of the small-amplitude breather, the system (|3.14[) 
is equivalent to the original system (|2.1[) . 

2. Despite the fact that the original 'PT-symmetric 
system includes gain and loss, the amplitude system was 
shown to be conservative. 

3. The breather solution was constructed as the 
asymptotic expansion (|4.12|t . (|4.13|> . (|4.14p . 

4. We have proved that all small-amplitude breathers 
are stable on the timescale t < e~ 2 . The small-amplitude 
breather was shown to be a "simple" oscillation — it 
cannot have any modulating frequencies in its spectrum. 

5. Breathers were shown to be common occurrences 
in the 'PT-symmetric chains of dimers. In particular, 
breathers are born in collisions of the low- and high- 
frequency solitons. 

In conclusion, we need to make three remarks. The 
first one is on the VT breathers versus conservative 
breathers and limit cycles. 

The 'PT-symmetric breathers are different from their 
conservative counterparts in that their associated physi- 
cal observables (e.g. energy and momentum) are not sta- 
tionary but oscillate in time. From this point of view, the 
VT breathers are similar to the time-periodic solitons in 
dissipative systems (57l - l59j . However there is an impor- 
tant distinction between the latter two categories too. 
Namely, the dissipative solitons are limit cycles (in an 
infinite-dimensional phase space); their amplitudes and 
periods are determined uniquely by the parameters of 
the system. On the contrary, the VT breathers arise as 
members of two-parameter families, similar to periodic 
trajectories in Hamiltonian systems. 

The second remark is on the radiation from the 
breather. Using the singular perturbation expansion, the 
breather can be constructed to any order in e. All higher- 
order corrections A n , B n are expressible as powers of Aq, 
Bq and decay to zero as |x| — > 00. There is no radiation 
to any order e™, n = 0, 1, 2, 3, .... 

However our simulations do reveal radiation waves 
from the breathers, with the amplitude of waves grow- 
ing as the amplitude of the breather is increased. The 
reason why the asymptotic expansion does not capture 
these waves is that the amplitude of radiation is expo- 
nentially small in e. (The exponential smallness does not 
imply that the radiation is invisible for finitely small e 
though.) 

The frequency of the radiation can be determined on 
the basis of standard considerations. Indeed, the spec- 
trum of linear excitations of the system (13. 2|) consists of 
two branches, u = k 2 + loq and u = k 2 — ujq [Fig. |BJ, 
while the breather of the amplitude e 1 / 2 has two basic 
frequencies, cj — e and — uj — e [see Eq. (|4.13|) ]. The 
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FIG. 6. The dispersion curves of the "PT-symmetric system 
(|2.1[) . The black dots indicate the two frequencies of the 
breather. The dashed line marks the frequency of radiation. 



term a 2 b* in (|3.2p oscillates at a combination frequency 
3wo— e which falls in the linear spectrum. Hence the dom- 
inant frequency of the resonant radiation will be 3wo — e, 
as indicated by the dashed line in Fig. [6l (Note that 
the frequency wo — e does not resonate with the bottom 
branch since the a and b modes are not coupled to the 
linear order.) 

Finally, we note that the breathers realise the periodic 
light switching between the waveguides with gain and 
loss. Unlike oscillations in structureless linear [IHIIil and 
nonlinear [2(| Hl| VT couplers, the breathers describe 
switching between spatially extended waveguides. Here, 
the nonlincarity suppresses the beam diffraction while 
the spatial coupling synchronizes the power oscillations 
across the beam. 



By rescaling P, Q, X and redefining e, we can always 
arrange that fj, = 1 in equations (14.41) : 

P" -P + 2(\P\ 2 + 2\Q\ 2 )P + — (|Q| 4 - 2\PQ\ 2 )P = 0, 

W 

Q" -vQ + 2{\Q\ 2 + 2\P\ 2 )Q + —{2\PQ\ 2 - \P\ A )Q = 0. 

Wo 

(Al) 

Note that we are not setting v equal to 1, along with p. 

For e = 0, the system (jAll) has even and odd solutions 
with n humps (n = 1,2,...), with both P and Q being 
nonzero |48|, l52j . Each of these can be used as a starting 
point in the regular perturbation expansion in powers of 
e. 

In particular, the solution of the system (IA1|) with e = 
0, with an even single-humped P(X) and an even single- 
humped Q(X), exists for a~ 2 < v < a 2 , where a 2 = 
|(a/I7 - l) 2 « 2.438, a- 2 « 0.410 0,^. Therefore 
the system (|Aip with sufficiently small nonzero e will 
also have a localised solution for any v between a~ 2 and 
a 2 . 

The solution with a two-humped even P(X) and a two- 
humped odd Q(X) exists for f3 2 < v < 1, where (3 2 = 

^(v^-a) 2 ^ o.3i5 Hi BH. 

All these soliton-like solutions of the system f|4.4[) give 
rise to breather solutions of the 'PT'-symmetric system 
(|2.1[) . Thus for each n > 1, the system (|2.1j) has a two- 
parameter family of nonpropagating breather solutions 
with n humps. Representatives of the family are different 
in the amplitude and width of the humps, as well as the 
contrast of the \u\ 2 - and \v\ 2 -oscillations. 



Appendix B: No internal modes for the 
small-amplitude breather 
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Appendix A: More general breather solutions 

In this Appendix, we briefly comment on other solu- 
tions of the system (|4.4p — more general than a nearly 
symmetric configuration (|4.7[) . (|4.6I) . (|4.10[) . 



The aim of this Appendix is to show that the operator 
(|5.13|) does not have any discrete eigenvalues. To this 
end, we note that if A =/= 0, the bottom component of 
(|5T31) gives 

J p 2 (X)zW(X)dX = 0, (Bl) 

where z^ ) = sechV is the null eigenvector of the operator 
L . The constraint (|B1[) defines a subspace of the space of 
square integrable functions; we will denote this subspace 
6. 

On the subspace S, the operator L$ is positive definite; 
hence we can write (|5.13[) as another scalar eigenvalue 
problem, alternative to (|5.15l) : 

L +P2 = -X 2 L^p 2: p 2 g 6. (B2) 

Assume the nonsymmetric matrix-differential operator 
in (I5.13[) has nonzero eigenvalues Ai,A2, The corre- 
sponding eigenvalues —A 2 < — A| < •■■ of (|B2[) are real, 



13 



and the associated eigenfunctions p 2 {X) can also be cho- 
sen real. The lowest eigenvalue can be found as the min- 
imum of the Rayleigh quotient: 



A-l = min 



(p 2 ,L + p 2 ) 
{P2,Lq 1 P 2 ) 



(B3) 



Since both L + and Lq are positive definite, Eq. (IB3|) im- 
plies that the eigenvalue — A^ of the generalised eigen- 
value problem (|B2[) is positive. Hence Ai lies in the 
gap of the continuous spectrum of the operator (|5. : 
Ai = Uji, — 1 < uji < 1. 

On the other hand, any function from © can be ex- 
panded over the continuous spectrum eigenfunctions of 
the operator Lq: 



P2 (X) = / Tl(k)z k (X)dk, 



(B4) 



where L Q Zk — (1 + 2k 2 )zk, — oo < k < oo. Writing 
L + as L + | sech 2 Af and substituting (|B4|) in (|B3[) . the 
Rayleigh quotient becomes 



/ TZ 2 {k)(l + 2k 2 )dk + jjpl sech 2 XdX 
fK 2 {k)(l + 2k 2 )~ 1 dk ' 



(B5) 



The first term in the numerator of (IB5[) is greater than 
the denominator; hence the quotient is greater than 1. 
This contradicts the fact that the eigenvalue Ai is in the 
gap of the continuous spectrum of the operator (|5.13|) . 
The contradiction proves that the operator (|5 . 13[) cannot 
have discrete eigenvalues. 
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